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SECTION  I 


INTRODUCTION 


In  netted  air  defense  and  air  traffic  control  systems,  data  from 
the  long  range  radars  are  routed  to  a  Sector  Operations  Center  and 
stereographically  projected  onto  a  common  coordinate  plane  for  presenta¬ 
tion  to  system  operators  on  the  display  consoles.  In  this  manner,  the 
overlap  coverage  of  the  radars  is  exploited  and  a  composite  air  sur¬ 
veillance  picture  is  presented. 

The  projection  of  aircraft  information  from  system  radars  may  be 
separated  into  two  parts  namely  conversion  and  transloruarloa  for 
convenience  in  analysis.  The  two  aspects  of  the  problem  may  be 
defined  as  follows; 

(a)  Conversion  of  sl.'Ut  range,  azimuth  and  height  data 
into  rectangular  coordinates  in  a  plane  with  the 
radar  site  as  origin. 

(b)  Transformation  of  these  coordinates  from  the  radar 
plane  into  rectangular  coordinates  in  the  common 
coordinate  plane. 

Reference  3  reports  on  an  accurate  and  simple  solution  to  the 
conversion  problem.  The  data  of  slant  range,  height  and  azimuth 
angle  are  stereographically  converted  to  rectangular  coordinates  in 
a  plane  tangent  to  the  earth  at  the  radar  site  location.  Assuming 
tangent  planes  have  been  so  established  at  multiple  radar  sites,  it 
is  desired  to  transform  the  coordinates  of  radar  data  in  the  various 
radar  planes  into  rectangular  coordinates  in  the  common  coordinate 
plane  with  accuracy  and  simplicity.  All  system  computations  such  as 
aircraft  tracking,  Interceptor  control  and  information  display  are 
performed  in  the  common  coordinate  plane. 

This  report  focuses  on  the  transformation  aspects  of  the  pro¬ 
jection  problem.  Its  purpose  is  to  supplement  Reference  3  and  to 
derive  and  develop  the  stereographic  transformation  equations  more 
completely  than  is  accomplished  in  References  1  and  2.  In 
particular  this  report: 

(a)  Derives  the  stereographic  transformation  equation 
in  complex  notation. 
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(b)  Derives  the  closed  form  solution  to  the  stereo¬ 
graphic  transformation  equation  in  terms  that 
are  suitable  for  real  time  computation. 

(c)  *Derive3  approximations  tc  the  closed  form  solu¬ 

tion  of  the  stereographic  transformation  equation 
in  terms  that  are  suitable  for  real  time  computation. 

(d)  Depicts  error  curves  assor  ,ated  with  the  approxima¬ 
tions  to  the  closed  form  solution  of  the  stereographic 
transformation  equation. 

(e)  Depicts  contours  of  constant  distortion  induced  by 
the  mapping  of  points  on  the  surface  of  the  earth 
ellipsoid  onto  the  surface  of  conformal  spheres  of 
varying  radii. 

(f)  Depicts  contours  of  constant  distortion  induced 
by  the  projection  of  the  surface  of  a  conformal 
sphere  onto  a  common  coordinate  plane. 

(g)  Presents  contours  of  constant  distortion  induced 
by  the  overall  projection  of  points  on  the  sur¬ 
face  of  the  earth  ellipsoid  onto  a  comon  coordi¬ 
nate  plane. 

(h)  Demonstrates  how  the  contours  of  constant  distor¬ 
tion  are  influenced  by  the  choice  of  an  earth's 
radius  and  reconmends  a  method  for  determining 
an  optimum  earth  radius  that  minimizes  the  dis¬ 
tortion  associated  with  the  overall  projection 
process. 

(i)  Demonstrates  the  Impact  of  tne  choice  of  the  origin 
of  the  common  coordinate  system  on  the  distortion 
contours. 

(j)  Summarizes  the  stereographic  cc aversion  and  trans¬ 
formation  equations. 


SECTION  II 

THE  TRANS FORMAT ION  EQUATION 


INTRODUCTION 

Appendix  A  derives  the  equation  for  transforming  the  coordinates 
of  data  from  the  plane  of  the  radar  site  to  a  common  coordinate  plane 
in  complex  notation.  This  section  of  the  report  presents  the  trans¬ 
formation  equation  in  complex  notation,  explains  the  various  terms 
therein  and  expresses  equation  (1)  in  terms  of  its  real  and  imaginary 
components.  Because  of  the  computational  difficulties  associated 
with  calculating  aircraft  locations  In  the  common  coordinate  plane 
via  the  derived  expression,  several  approximations  that  combine  sim¬ 
plicity  with  accuracy  are  also  presented. 

From  equation  (A-13)  of  Appendix  A,  the  transformation  equation 
in  complex  notation  is: 


W 


Ze~16  +  W 


ZV  e-1B 
r 


1  - 


4E 


(1) 


where: 

W  -  U  +  IV  are  the  rectangular  coordinates  of  an 
aircraft  with  respect  to  the  common  coordinate 
origin. 

Z  ••  X  +  iY  are  the  rectangular  coordinates  of  an 
aircraft  with  respect  to  the  radar  coordinate 
origin. 

3  ■  is  a  rotation  angle  which  ackes  the  axes  of 

the  Z  and  W,  the  radar  and  common  coordinate 
planes  more  nearly  parallel. 

■  Ur  +  1V_  are  the  rectangular  coordinates  of  a 
radar  with  respect  to  the  common  coordinate 
origin. 
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■  0^  -  iVr  is  the  complex  conjugate  of  Wf. 

E  -  the  radius  of  a  spherical  earth. 

With  respect  to  the  components  of  the  W  term  in  equation  (1) , 
the  positive  U  axis  is  directed  toward  east  and  the  positive  V  axis 
is  directed  toward  north  at  the  origin  of  the  system.  With  respect 
to  the  components  of  the  Z  term  in  equation  (1) ,  the  positive  X  axis 
is  directed  toward  east  and  the  positive  Y  axis  is  directed  toward 
north  at  the  origin  of  the  radar  coordinate  system.  The  coordinate 
axes  in  the  radar  coordinate  plane  are  shown  in  Figure  1. 

Given  the  measured  slant  range  from  the  radar  to  an  aircraft , 
the  height  of  an  aircraft  above  sea  level  and  the  elevation  of  the 
radar  above  sea  level,  the  stereographic  ground  range  R,  the  dis¬ 
tance  from  the  radar  to  the  aircraft  in  the  radar  plane,  is  obtained 
as  described  in  Reference  3.  Having  obtained  the  ground  range  R, 
the  rectangular  coordinates  of  the  aircraft  with  respect  to  the 
radar  coordinate  origin  may  be  expressed  as: 

Z  »  |z|ei^ir/2)"01  -  KfiiI(ir/2)-e]  M  R[Sln  9  +  ±  Cos  ej 


Therefore,  the  rectangular  coordinates  of  Z  may  be  written  as: 

X  -=  R  Sin  0 
Y  -  R  Cos  0 


where : 

R  is  the  projected  ground  range  of  the  aircraft. 

6  is  the  azimuth  angle  measured  clockwise  from  the 
positive  Y  axis  which  is  oriented  toward  north 
at  the  location  of  the  radar  site. 

The  ang3e  0  in  equation  (1)  is  derived  in  equation  (A-9)  of 
Appendix  A  and  depends  simply  upon  the  positional  coordinates  of 
the  radar  site  with  respect  to  the  origin  of  the  common  coordinate 
plane.  These  coordinates  are  defined  by  their  latitudes  L  and 


longitudes  X.  Latitudes  north  of  the  equator  are  defined  as  positive 
and  latitudes  south  of  the  equator  are  defined  as  negative.  Longi¬ 
tudes  are  defined  as  positive  east  of  the  prime  meridian  and  negative 
west  of  the  prime  meridian.  From  equation  (A-9) ,  S,  which  is  a  con¬ 
stant  for  any  radar  with  respect  to  a  particular  conmon  coordinate 
origin,  is  defined  as  follows: 


Tan 


-1 


(Sin  L  +  Sin  L  )  Sin  (X  -X) 
o  o 


Cos  L  Cos  L  +  (1  +  Sin  L  Sin  L  )  Cos  (X-X  ) 
o'  o  o 


(2) 


where: 


L,X  are  the  latitude  and  longitude  of  the  origin 
of  the  tangent  plane  at  the  radar  site. 


L  ,X 
o  o 


are  the  latitude  and  longitude  of  the  origin 
of  the  conmon  coordinate  plane. 


The  angle  6  may  be  considered  as  a  rotation  of  the  Z  plane  with 
respect  to  the  W  plane.  The  effect  of  the  rotation,  which  ia 
counterclockwise  when  X  <  XQ  and  clockwise  when  X  >  X0,  is  to  make 
the  axes  of  the  radar  plane  more  nearly  parallel  to  the  axes  of 
the  conmon  coordinate  plane. 

Considering  which  represents  the  coordinates  of  a  radar  site 
with  respect  to  the  common  coordinate  origin.  Appendix  B  derives  the 
equations  which  stereographically  map  points  on  the  surface  of  the 
earth  onto  a  rectangular  coordinate  plane  for  a  spherical  earth. 

These  equations  are  used  to  project  points  such  as  the  location  of 
a  radar  site  onto  the  common  coordinate  plane.  From  equation  (B-5) 
of  Appendix  B,  the  rectangular  coordinates  of  a  radar  site  located 
at  latitude  L  and  longitude  X  with  respect  to  the  common  coordinate 
plane  are: 


2E 


2E 


Sin  (X-XQ)  Cos  L 


1  +  Sin  L  Sin  1^  +  Cos  L  Cos  1^  Cos  (X-Xo> 

Sin  L  Cos  -  Cos  L  Sin  Cos  (X-X^ 

1  +  Sin  L  Sin  1^,  +  Cos  L  Cos  lQ  Cos  (X-XQ) 


(3) 
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where  * 

L,X  are  the  latitude  and  longitude  of  the  radar 
site,  the  coordinates  of  which  are  projected 
onto  the  common  coordinate  plane. 

Lo,Xq  are  the  latitude  and  longitude  of  the  origin 
of  the  coeuon  coordinate  plane. 

are  the  rectangular  coordinates  of  the  point 
designated  by  L,X  in  the  common  coordinate 
plane. 

Thus  far,  a  brief  description  of  each  of  the  terms  in  equation  (1) 
has  been  presented.  Equations  for  obtaining  the  rectangular  coordi¬ 
nates  of  an  aircraft  in  a  plane  tangent  to  the  earth  at  the  location 
of  a  radar  site  have  been  presented.  Assuming  a  common  coordinate 
tangent  plane  has  been  established  at  a  location  different  from  the 
radar  site,  results  have  been  derived  for  the  rotation  angle  3  which 
make  the  axes  of  the  radar  plane  more  nearly  parallel  to  those  of 
the  comnon  coordinate  plane.  Additionally,  equations  have  been 
derived  which  project  the  earth  coordinates  of  a  radar  onto  a  point 
in  the  common  coordinate  plane.  It  is  now  desired  to  employ  this 
information  in  expressing  equation  (1)  in  terms  of  its  real  and 
imaginary  components. 


CLOSED  FORM  TRANSFORMATION  EQUATION 

Equation  (1)  is  the  complex  representation  of  an  aircraft's  posi¬ 
tion  in  the  common  coordinate  plane.  Expressing  the  numerator  of  this 
equation  in  terms  of  R,  0,  6  and  |wr|,  we  note  that: 


Ze~ie  -  (X  +  iY) (Cos  3  -iSin  3)  -  (X  Cos  3  +  Y  Sin  8)  +  i(Y  Cos  8-X  Sin  $) 
-  R  Sin  (6  +  3)  +  iR  Cos  (0  +  6)  *  A  +  IB 

Therefore,  the  numeratev  of  equation  (1)  is: 

[R  Sin  <8  +  8)  +  |Wr|  Sin  y]  +  i[R  Cos  <9  +  3)  +  |wj  Cos  y]  -  C  +  ID 
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The  denominator  of  equation  (1)  may  be  expressed  in  terms  of  R, 
0,  $,  and  E  as  follows  * 


Employing  the  above  definitions,  equation  (1)  may  be  written  as: 


W  -  U  +  iV  -  C  +  10  «  (GC  -  FD)  +  i(GD  +  FC) 

G  -  iF  g2  +  F2 


with: 


GC  ~  FD 

2  2 
G  +  F* 


and 


GD  4-  FC 

2  2 
G*  +  F 


2  2 

Evaluating  GC  -  FD,  GD  +  FC  ,  and  G  +  F  we  obtain: 


GC  -  FD  -  R  Sin  (6  +  8)  +  U  - 

r 


UrR2  +  |wJ2R  Sin  [2y  -  (6  +  8)] 
4E2 


GD  +  FC  -  R  Coa  (0  +  0)  +  V  - 

r 


Vr  +  lWr'  R  CoB  [2y  ~  (e  s'  B)] 
4E2 

2„2 


2  2  2 1 W  | R  Cos  [y  —  (9  +  $)]  |wJzr' 

3  +  f  -  i  — - - = - +  ...y ,  , 


4E 


(4E2)2 


Therefore,  the  closed  form  solution  to  the  stereograohic  trans¬ 
formation  equation  is: 


R  Sin  (6  +  8)  +  Ur  - 


U  - 


R2  +  | Wr | 2  R  Sin  [2y  -  (9  +  B)] 
4£2 


1  - 


2  wjR  Cob  ty  -  (0  +  8)] 

'|Wrl  R 

2 

2  + 
4EZ 

•  <4E2)  . 

R  Cos  (9  +  8)  +  Vf  - 


V  - 


Vr  R2  +  |Wr | 2  R  Cos  [2y  -  (9  +  8)] 
4E2 


1  - 


2 J W  J R  Cos  [y  -  <8  +  8)) 


4E 


W  R 
r 


-.2 


4E 


(4) 
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Because  of  the  computational  difficulties  associated  with  cal¬ 
culating  aircraft  locations  via  these  expressions  for  each  radar 
datum,  an  approximation  that  combines  suitable  accuracy  with  minimal 
processing  requirements  is  desirable. 


APPROXIMATE  TRANSFORMATION  EQUATION 

Equation  (1),  which  for  convenience  is  reiterated  below,  projects 
aircraft  onto  a  tangent  plane  centered  at  the  common  coordinate  origin. 


W 


where; 


Z  +  W 


z  w 


1  - 


4E 


Z  -  Ze"10 


(5) 


Expanding  equation  (5)  yields: 

,  ,-1 


W  -  (Z  +  wr) 


Z  W 


1  - 


4E 


(Z  +  wr) 


1  + 


Sftr) 


w  +  w 


r  r 


n*l 


• _  1  * 

n 

00 

-  —  f  - 

w  z 

r 

t  « 

+  Z  +  Z 

y; 

W  Z 
r 

4E2 

•  - 

n-l 

4E2 

"  (w_)n_1  (zV 


(W  )n  (z')1^1 


,  yc  )  ,  y»  /  y&  ) 

wgf|wg|2L-r  2n —  +z  +  — rs 

r  r  n-l  (4E2)n  n-l  (4E2)n 


(6) 


U  +  - 5-  Zw  - 


2  "  (wr)n  (z')n+1  ^  (vr)n  (z’)n+1 


r  4E  n-0  (4EZ) 


2.n 


n-0  (4EZ)n 


W  + 
r 


l»rl2' 

:v<v”  (z')0+1  ’ 

1  +  2 
4E* 

n-0  (4E2)” 

Let  d  and  y  be  the  azlnuths  of  Z  and  Wr  respectively.  Positive 
values  of  6  are  measured  clockwise  from  true  north  at  the  radar  site 
and  in  the  radar  plane.  Positive  values  of  y  are  measured  clockwise 
in  the  common  coordinate  plane  relative  to  true  north  at  the  origin 
of  the  tf  plane.  Positive  values  of  the  rotation  angle  B,  which  is 
derived  in  Appendix  A  and  presented  in  equation  (2),  are  measured 
counterclockwise  from  east  at  the  radar  site  and  in  the  radar  plane. 
In  order  that  9,  y  and  3  may  properly  be  combined ,  it  is  necessary 
to  describe  6  and  y  in  the  same  coordinate  system  as  that  in  which  B 
is  measured.  The  angles  corresponding  to  6  and  y  in  a  system  where 
positive  angles  are  measured  counterclockwise  from  east  are  t(*/2)-9] 
and  [ (tr/2)  — y J  respectively.  Therefore,  we  may  write: 


(H.)"  -  |w  |"  -  TTl 


(z')"+1  -  (Ze-16) 


n+1 


iz|eif<*/2>  -  e]e-iB 


n+1 


where: 


|  ^  |  n+1  ei(n+l)  [(ir/2)  -  (0  +  B)]  _  i  z  i n+X  ei(n+l)  [  (ir/2>  -  6  J 


e  -  (e  +  p) 


From  equation  (6)  and  employing  the  above  definitions,  U  may  be 
written  as: 


W  -  W  +  K 
r 


Wr  +  K 


in  -|  n+1  i 

V  Z|  -in[(w/2)  -  y]  _i(n+X)  [(w/2)  -  0  ] 

- e  e 


iTo  <4E2)n 


Oft  | 


in  , 2 , n+1 


n-0  (4E2)n 


i[ (tt/2)  -  {o'  +  0(9*  -  y) >1 


Wr  +  K 


y  |wt[°  Izl"*1 

n-0  <4E2>° 


Sin  f  (9+B)+n (9+6-y)  ]  +  i  Cos  t  <6+B>++i (6-+6-y )  ] 
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Therefore,  the  components  of  U  are: 


V 


U  +  K 
r 


y-  [wr|°|Z|ntl 

n-0  <4E2)n 


Sin  1(6  +  g)  +  n  (0  +  0  -  y)1 


V 


V  lwri  izl 

V  +  K  2—i  - T- -  Cos  [(6  +  0)  +  n  <9  +  8  -  y>] 

n-0  <4EV 


where: 


K  - 


1  + 


Approximations  to  U  and  V  may  be  obtained  by  terminating  this  series 
solution  to  equation  (1)  after  an  appropriate  number  of  terms. 


First  Order  Approximation 

The  first  order  approximation  to  equation  (1)  is  obtained  by 
taking  the  first  term  of  equation  (7)  as  follows: 

U  -  Uf  +  K  | Z |  Sin  (0  +  B)  -  Ur  +  K  R  Sin  (0  +  B) 

V  -  V  +  K  |  Z I  Cos  (6  +  B)  -  V  +  K  R  Cos  (6  +  B) 

r  r 

The  error  in  this  approximation  to  equation  (1)  is  on  the  order 
of  the  first  term  neglected.  Values  for  the  first  term  neglected, 
namely  k|w_Ir2/4E  ,  are  nresented  in  Figure  2  for  ground  ranges  from 
25  to  200  miles  and  for  ]Wr|  from  0  to  550  miles. 

Figure  2  shows  the  transformation  errors  associated  with  a  first 
order  approximation  to  equation  (1).  These  errors  exceed  0.4  miles  at 
large  values  of  R  and  | 1  end  approach  0.25  miles  at  intermediate 

values  of  R  and  |Wr|.  The  composite  projection  error  is  the  sum  of 

the  transformation  error  as  shown  in  Figure  2  and  the  conversion  error 
as  shown  in  Reference  3. 


When  aircraft  height  Information  is  available,  the  upper  limit 
on  the  conversion  error  is  expected  to  be  0.1  miles.  This  result  has 
been  obtained  by  consideration  of  Figure  6  of  Reference  3.  Assuming 
a  maximum  distance  between  a  radar  site  and  the  projection  center  of 
275  miles,  the  maximum  transformation  error  from  Figure  2  herein  is 
anticipated  to  be  0.25  miles  using  a  first  order  approximation.  The 
composite  maximum  error  is  0.35  miles,  which  is  deemed  unacceptable. 
Consequently,  when  aircraft  height  information  is  available,  a  higher 
order  approximation  to  equation  (1)  that  will  essentially  decrease 
the  transformation  error  to  zero  as  will  be  shown  hereafter,  is 
required . 

From  Figure  12  of  Reference  3,  it  is  noted  that  errors  in  excess 
of  one  mile  occur  in  converting  slant  range  to  ground  range  when  air¬ 
craft  height  data  is  unavailable.  The  fundamental  remedy  to  this 
problem  is  to  obtain  height  information  on  the  aircraft  and  thereby 
diminish  the  relatively  large  conversion  error  rather  than  attempting 
to  decrease  a  first  order  transformation  error  that  may  at  most  equal 
0.25  miles.  Therefore,  when  aircraft  height  Information  is  unavail¬ 
able,  a  first  order  approximation  to  equation  (1)  is  acceptable. 


Second  Order  Approximation 

The  second  order  approximation  to  equation  (1)  is  obtained  by 
taking  the  first  two  terms  of  equation  (7)  as  follows: 


U  -  U  +  K 
r 


R  |W |  Sin  [2(6  +3)  -  y] 
R  Sin  (8  +  6)  +  - - - s - 


4E 


V  *  V  +  K 
r 


R  Cos  <6  +  3)  + 


R  | Wr !  Cos  [2(6  +  3)  -  Y) 


4  E 


The  error  in  this  approximation  is  on  the  order  of  the  third 
term  of  the  infinite  series.  Values  for  the  third  term,  namely 
K|wr |2R3/16E*  are  presented  in  Figure  3  for  ground  ranges  from  75  to 
200  miles  and  for  JW_|  from  0  to  550  miles.  Examination  of  Figure  3 
shows  that  the  second  order  approximation  to  equation  (1)  essentially 
yields  identical  results  to  the  closed  form  solution  to  the  stereo- 
graphic  transformation  equation. 


For  programing  convenience,  equations  (9)  may  be  further  sim¬ 
plified  by  a  numerical  analysis  of  the  second  order  range  term.  When 
height  data  is  available  a  suitable  expression  for  the  ground  range  R, 
which  is  derived  in  equation  (8)  of  Reference  3,  is  given  by: 


R 


[S2  -  (H-h)2ft 

1  +  f 


where : 

S  is  the  measured  slant  range  of  the  aircraft. 

H  1e  the  elevation  of  the  aircraft  above  seu  level. 

h  is  the  elevation  of  the  radar  above  sea  level. 

E  is  the  radius  of  a  spherical  earth. 

Substituting  this  result  into  the  U  component  portion  of  equation  (9) 
yields : 


U 


U  +  K 

c 


(S7-(H-h)2ft 


Sin  (6+6)  + 


[S2-(H-h)21  [WrlSi..  [2(6+6)  -  y] 

2  2  2 
(1+f)2  4E2 


Calculating  the  maximum  magnitude  of  the  {|w„  Sin  [2(6+6)  -  y]}/4E2) 
term,  we  obtain  the  result  1.056  x  10“^  when'  Wr|  equals  500  miles. 
This  multiplier  of  the  Rz  term  is  miniscule;  consequently,  the 
denominator  of  the  R2  term  may  be  changed  to  (1  +  2/E)  while  retain¬ 
ing  an  accuracy  of  better  than  3.1  x  10“^  miles  for  values  of  R 
through  225  miles  and  for  values  of  |Wr|  through  500  miles. 


Rewriting  equations  (9)  as  suggested,  we  obtain  a  revised  approx¬ 
imation  to  equation  (1)  that  retains  the  accuracy  of  equations  (9) 
while  simplifying  the  formulation  and  processing  requirements.  The 
simplified  second  order  approximation  is: 


U 


V 


U.  +  T 


Vr  +  T 


tS^Ol-h)2]1*  Sin 


JS^Oi-h)2]5*  Cos 


(6+6)  +  tS2-(H-h)2]  l»rl  Sin  [2(6+6)  -  y1 

4E2 

(6+6)  +  [S2-(H-h)2l  |Wr|  Cos  [2(6+6)  -  y] 

4E2 


(9a) 
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where: 


1  + 


4E2 


1+l 


T  is  a  site  dependent  constant  and  say  be  precalculated  and  stored  as 
an  adaptation  parameter . 
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SECTION  III 


DISTORTION  INDUCED  BY  PROJECTION 


INTRODUCTION 

The  earth  is  not  a  sphere  but  is  an  ellipsoid  that  is  approxi¬ 
mately  23  miles  greater  in  equatorial  than  polar  diameter.  This  fact 
has  been  deliberately  ignored  until  this  section  of  the  report  as  it 
appeared  wiser  to  Initially  consider  a  spherical  earth  and  later  to 
consider  the  implications  of  the  earth's  ellipticity.  A  topic  that 
must  now  be  addressed  is  the  value  of  an  earth's  radius  that  should 
be  used  in  any  of  the  equations  which  are  dependent  upon  the  earth's 
size  and  shape.  A  related  topic  is  the  distortion  or  error  that  is 
introduced  by  the  entire  mapping  process. 

Any  projection  of  an  ellipsoid  onto  a  plane  distorts  either  dis¬ 
tances,  angles  or  areas,  it  being  impossible  to  simultaneously  pre¬ 
serve  all  of  these  parameters.  Furthermore,  no  projection  of  an 
ellipsoid  onto  a  plane  can  simultaneously  preserve  both  angles  and 
areas.  The  stereographic  projection  preserves  angles  but  distorts 
distances. 

This  section  of  the  report  will  address  the  manner  in  which  the 
earth's  ellipticity  must  be  taken  into  account.  Additionally,  the 
magnitude  of  the  distortion  induced  by  the  projection  of  an  ellipsoid 
onto  a  plane  will  be  discussed.  Considerations  relative  to  the  choice 
of  an  earth's  radius  in  the  transformation  equations  will  be  discussed. 
The  effect  of  the  choice  of  an  earth's  radius  upon  the  overall  distor¬ 
tion  associated  with  stereographic  projection  will  be  quantitatively 
presented. 

By  way  of  introduction  to  these  topics,  the  stereographic  pro¬ 
jection  must  be  considered  in  greater  detail. 

STEREOGRAPHIC  PROJECTION 

The  process  of  stereographically  projecting  points  onto  a  plane 
is  accomplished  in  what  may  be  considered  as  two  distinct  steps. 

The  surface  of  the  earth  is  approximated  by  the  surface  of  an  ellip¬ 
soid,  which  is  obtained  by  rotating  an  ellipse  about  its  minor  axis. 
The  first  step  maps  points  on  or  above  a  designated  point  on  the 
ellipsoid  onto  a  sphere  which  will  be  called  the  conformal  sphere. 

The  8  cond  step  maps  points  on  the  conformal  sphere  onto  a  plane 
tangent  to  the  conformal  sphere. 
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The  mapping  of  points  on  or  above  a  designated  location  on  the 
ellipsoid  onto  a  sphere  is  accomplished  by  the  following  relations 
as  indicated  on  page  34  of  Reference  4: 


-  e  Sin  L 
+  e  Sin  L 


where: 

L, A  are  the  latitude  and  longitude  of  the  point  on  the 
ellipsoid. 

t 

4, A  are  the  latitude  and  longitude  of  the  corresponding 
point  projected  onto  the  sphere.  This  latitude  is 
defined  as  the  conformal  latitude. 


e  *  .08199189  is  the  eccentricity  of  the  earth. 


Therefore,  the  equations  for  mapping  an  ellipsoid  onto  a  conformal 
sphere  and  for  stereographically  mapping  the  conformal  sphere  onto  a 
plane  are  the  same  as  those  that  have  been  presented  for  a  spherical 
earth  except  that  the  conformal  latitude  $  must  be  used  in  place  of 
the  geographic  latitude  L.  For  computational  convenience,  the  con¬ 
formal  latitude  $  may  be  obtained  from  the  geographic  latitude  L  by 
a  series  approximation  to  equation  (10) . 

As  noted  on  page  34  of  Reference  4,  the  mapping  from  the  ellip¬ 
soid  to  the  sphere  is  conformal.  The  mapping  from  the  sphere  to  the 
plane  is  also  conformal.  Each  step  preserves  angles;  hence,  the 
resultant  stereographic  projection  is  a  conformal  representation  of 
the  earth  ellipsoid  on  a  plane. 


(10) 


DISTORTION  AND  EARTH'S  RADIUS 

The  surface  of  an  ellipsoid  cannot  be  molded  into  a  sphere  with¬ 
out  some  distortion  or  error.  Similarly,  the  surface  of  an  appropri¬ 
ately  chosen  sphere  cannot  be  spread  out  flat  without  some  distortion. 

In  general,  arc  lengths  that  are  projected  onto  a  planar  surface  are 
stretched  as  the  distance  between  the  point  of  tangency  and  the  point 
being  projected  increases.  The  magnitude  of  the  stretching  is  non¬ 
linear.  Therefore,  groups  of  positional  reports  from  the  same  air¬ 
craft  exhibit  distortions  when  projected  on  a  plane.  Data  from  a  con¬ 
stant  speed  aircraft,  that  is  constrained  to  move  at  a  constant  altitude 
above  the  earth's  surface,  portray  apparent  accelerations  when  projected 
onto  a  plane. 
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Intuitively,  it  appears  evident  that  the  distortions  described 
above  can  be  minimized  by  a  judicious  choice  of  the  earth's  radius. 

The  following  paragraphs  will  show  quantitative  relationships  between 
the  choice  of  an  earth's  radius  and  the  resultant  distortion  for  the 
mapping  from  the  ellipsoid  to  the  conformal  sphere  and  for  the  pro¬ 
jection  of  the  conformal  sphere  onto  a  plane. 

Mapping  from  Ellipsoid  to  Sphere 

The  scale  factor  associated  with  the  mapping  of  points  from  an 
ellipsoid  onto  a  sphere  is  defined  as  the  ratio  of  the  arc  length 
along  the  surface  of  the  sphere  to  the  corresponding  arc  length  along 
the  surface  of  the  ellipsoid.  A  scale  factor  of  unity  indicates  that 
distances  are  completely  preserved  and  that  there  is  no  distortion  in 
the  mapping  process.  From  page  86  of  Reference  5,  the  scale  factor 
k. ,  associated  with  the  projection  of  the  ellipsoid  onto  the  conformal 
sphere  is: 

,  _  ffl  _  E  Cos  j 
*1  "  ds2  “  N  Cos  L 

where : 

is  the  scale  factor. 

ds^  is  the  arc  length  along  the  surface  of  the  conformal 
sphere . 

d&2  is  the  arc  length  along  the  surface  of  the  ellipsoid. 

E  is  the  radius  of  the  conformal  sphere. 

L  is  the  geographic  latitude.  This  is  the  angle  which 
a  normal  to  the  surface  of  an  ellipsoid  makes  with 
the  equatorial  plane  as  shown  in  Figure  4. 

$  is  the  conformal  latitude  of  the  point  whose  geo¬ 
graphic  latitude  is  L. 

N  is  the  distance  HP  in  Figure  4  where  a  normal  to  the 

ellipsoid  at  point  M  has  been  drawn  until  it  Intersects 
the  minor  axis  at  P.  The  noraalN  is  mathematically 
expressed  as  Eq  (1  —  Sin^  L)"^  where  Eq  “  3444.054 
miles  is  the  earth's  equatorial  radius  and  where  - 
.00672267  is  the  earth's  eccentricity  squared. 


(ID 
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Figure  4 

NORMAL  TO  THE  ELLIPSOID 


Equation  (II)  expresses  the  scale  factor  k^  in  terns  of  the 
earth's  radius  at  the  equator,  the  earth's  eccentricity,  the  geo¬ 
graphic  latitude  and  the  radius  of  the  conformal  sphere.  Therefore, 
we  may  set  the  scale  factor  k^  equal  to  unity  at  a  designated  latitude 
and  solve  for  the  resultant  radius  of  the  conformal  sphere  according 
to  the  relationship: 


N  Cos  L 

E  »  — - —— 

Cos  $ 

Equation  (12)  determines  the  radius  of  the  conformal  sphere 
that  will  result  in  a  unity  scale  factor  at  the  latitude  that  is 
used  in  the  computation.  Table  1  shows  the  radii  of  the  conformal 
spheres  that  result  from  solving  equation  (12)  for  latitudes  L  of  0, 
20  and  60  degrees.  The  conformal  latitudes  <p  that  co-respond  to 
these  geographic  latitudes  are  also  shown. 


(12) 


21 


i 


Table  1 


Earth  Radius  Versus  Latitude  (Ellipsoid  To  Sphere) 


Geographic  Latitude 

Conformal  Latitude 

Earth  Radius 

0.0 

0.00 

3444.05 

20.0 

19.88 

3442.71 

40.0 

39.81 

3439.29 

Table  I  shows  the  radii  of  the  conformal  spheres  that  result 
from  solving  equation  (12)  using  the  indicated  geographic  latitudes. 
When  the  mapping  of  the  ellipsoid  onto  the  sphere  is  centered  at 
the  equator,  the  optimum  radius  of  the  conformal  sphere  is  the  earth's 
equatorial  radius  of  3444.05  miles.  When  the  mapping  Is  centered  at 
more  northerly  latitudes,  the  optimum  earth  radii  decrease  as  shown 
in  the  table. 

Let  us  now  determine  the  scale  factors  or  magnifications  that 
result  from  solving  equation  (11)  using  each  of  the  earth  radii  con¬ 
tained  in  Table  I  over  a  designated  set  cf  latitudes.  It  1b  evident 
that  the  radius  of  the  conformal  sphere,  which  is  calculated  based 
upon  a  particular  latitude,  will  result  in  magnifications  closest 
to  unity  over  a  geographical  area  that  1b  centered  at  the  same 
latitude.  To  quantify  this  judgment,  Table  II  shows  the  solutions 
to  equation  (11)  for  latitudes  between  36  and  44  degrees  when  E 
assumes  each  of  the  three  values  shown  in  Table  I. 

Table  II 

Magnification  Values 


Latitude/Corresponding  Radius  of  Conformal  Sphere 


Latitude 

0.0/3444.05 

20.0/3442.71 

40/3439.29 

36 

1.00116 

1.00077 

0.99977 

38 

1.00127 

1.00088 

0.99989 

40 

1.00138 

1.00099 

1.00000 

42 

1.00150 

1.00111 

1.00012 

44 

1.00162 

1.00123 

1.00023 
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Each  column  of  Table  II  shows  the  magnification  value  that  re- 
suits  from  solving  equation  (11)  using  the  indicated  value  of  the 
radius  of  the  conformal  sphere  while  varying  the  latitude  from  36 
to  44  degrees.  When  the  radius  of  the  conformal  sphere  is  calculated 
based  upon  unity  magnification  at  latitudes  of  either  0  or  20  degrees, 
distances  on  the  ellipsoid  over  all  latitudes  from  36  to  44  degrees 
are  stretched  when  projected  onto  the  respective  spheres.  When  the 
radius  of  the  conformal  sphere  is  based  upon  unity  magnification  at 
a  latitude  of  40  degrees,  distances  on  the.  ellipsoid  at  latitudes 
above  40  degrees  are  elongated  whereas  distances  below  a  latitude 
of  40  degrees  are  contracted  when  projected  onto  the  sphere.  As 
expected,  the  magnitudes  of  the  elongations  and  contractions  are 
minuscule.  Figure  5  depicts  the  situation  when  the  earth's  radius 
is  chosen  based  upon  the  latitudinal  center  of  the  geographical  area 
of  interest,  designated  L.  As  Indicated  in  Table  II,  this  radius 
of  the  conformal  sphere  is  optimum  in  that  it  results  in  magnifica¬ 
tions  closest  to  unity  over  the  area  of  interest. 

Mapping  from  Ellipsoid  to  Plane 

The  scale  factor  k_  that  is  associated  with  the  process  of  map¬ 
ping  points  on  the  conformal  sphere  onto  points  on  the  common  coordi¬ 
nate  plane  is  from  page  41  of  Reference  4: 


v  _  ,  - 

2  1  +  Sin  $  Sin  <J>  +  Cos  $  Cos  4>  Cos  (X-A  ) 

o  o  o 

where : 

<J>,A  are  the  conformal  latitude  and  longitude  of  the 
point  to  be  projected, 

$  ,A  are  the  conformal  latitude  and  longitude  of  the 
°  origin  of  the ’common  coordinate  plane.  This 

origin  will  be  referred  to  as  the  projection 
center , 

From  equation  (B-l)  of  Appendix  B,  it  will  be  noted  that  the  de¬ 
nominator  of  equation  (13)  may  be  vrltten  as: 

1  +  Sin  $  Sin  $  +  Cos  4  Cos  $  Cos  (A-A  )  ■  1  +  Cos  D 

O  0  0 

where : 

D  is  the  angle  that  subtends  the  great  circle  arc  from  ($,A) 

t°  <vv- 


(13) 
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The  magnification  ±2  equals  unity  at  the  projection  center  and 
increases  as  the  distance  from  the  projection  center  increases.  Con¬ 
tours  of  equal  scale  magnification  are  concentric  circles  about  the 
projection  center.  Figure  6  shows  contours  of  equal  magnification  for 
several  values  of  k2  where  the  horizontal  axis  is  a  range  of  longitudes 
between  -9  and  +3  degrees  and  where  the  vertical  axis  1b  a  range  of 
latitudes  between  36  and  44  degrees.  From  Figure  6,  it  will  be 
observed  that  the  scale  is  unity  at  the  projection  center,  that  the 
scale  increases  with  the  distance  from  the  projection  center  and  that 
contours  of  equivalent  magnification  are  symmetric  about  the  vertical 
axis.  The  projection  center  Lq ,  Xq  in  Figure  6  is  40°,  -3°. 

The  total  scale  factor  associated  with  the  projection  of  the 
ellipse  onto  the  plana  is  the  product  of  the  scale  factors  associated 
with  the  two  steps  of  the  mapping  process.  Therefore,  the  total 
magnification  designated  kfc  is: 

_ _ 2E  Cos  6 _ 

t  *  N  Cos  L  [1  +  Sin  $  Sin  4  +  Cos  4  Cos  4  Cos  (X-X  )j 

o  o  o 

Table  II  has  presented  the  magnification  values  associated  with 
the  projection  of  the  ellipsoid  onto  conformal  spheres  with  varying 
radii.  Figure  6  has  presented  contours  of  equal  magnification  for 
the  mapping  from  the  conformal  sphere  onto  the  plane.  Equation  (14) 
provides  a  method  for  combining  the  magnifications  associated  with 
the  two  ateps  of  the  mapping  process  over  a  designated  range  of 
latitudes  and  longitudes  for  a  specified  radius  of  the  conformal 
sphere . 

For  comparative  purposes  equation  (14)  has  been  evaluated  using 
as  the  radius  of  the  conformal  sphere: 

(a)  The  earth's  equatorial  radius 

(b)  The  earth  radius  that  is  optimized  for  the  mapping 
from  the  ellipsoid  to  the  sphere  when  the  mapping  is 
latitudinally  centered  at  40  degrees. 

Using  the  earth's  equatorial  radius  as  the  radius  of  the  conformal 
sphere,  Figure  7  shows  contours  of  equal  magnification  for  the  same 
range  of  latitudes  and  longitudes  that  were  used  in  Figure  6.  Examina¬ 
tion  of  Figure  7  shows  that  contours  of  equal  mangification  remain 
symmetric  about  the  vertical  axis.  Due  to  use  of  an  earth's  radius 
that  is  better  suited  to  a  projection  centered  at  the  equator  than  to 
a  mapping  centered  at  a  latitude  of  40  degrees,  the  contours  of  equal 
magnification  are  larger  than  is  necessary.  The  magnification  is 
1.00138  at  the  projection  center  and  lncraeees  away  from  the  projec¬ 
tion  center. 


(14) 
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Figure  6  CONFORMAL  SPHERE  TO  PLANE  MAGNIFICATION 
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Figure  8  shows  contours  of  equal  magnification  using  the  earth's 
radius  that  is  optimized  for  the  mapping  from  the  ellipsoid  to  the 
sphere  when  the  mapping  is  latltudinally  centered  at  40  degrees. 

With  cursory  consideration,  this  choice  of  an  earth's  radius  appears 
to  minimize  distance  distortions  over  the  geographical  area  of  interest 
since  is  unity  at  the  latitudinal  center  of  the  projection  and 
since  is  completely  determined  for  any  point  (L,X)  given  the  pro¬ 
jection  center  (L  ,  X  ). 

o  o 

The  magnification  contours  of  Figure  8  are  analogous  to  those 
of  Figure  6  except  that  the  contours  of  Figure  8  are  no  longer 
symmetric  about  the  horizontal  axis  through  the  projection  center; 
these  contours  are  flattened  to  the  north  and  elongated  to  the  south 
of  the  projection  center.  Like  Figure  6,  the  magnification  is  unity 
at  the  projection  center,  increases  away  from  the  center,  and  is 
symmetric  about  the  vertical  axis. 

Through  this  point  in  the  discussion.  Figure  6  has  presented 
contours  of  equal  magnification  associated  with  the  mapping  from  the 
conformal  sphere  to  the  plane.  These  contours  are  independent  of 
the  choice  of  an  earth  radius.  Figures  7  and  8  have  presented  con¬ 
tours  of  equal  magnification  based  upon  the  overall  projection  from 
the  ellipsoid  to  the  plane.  These  contours  are  dependent  upon  the 
choice  of  an  earth's  radius  and  have  been  depicted  for  radii  that 
minimize  the  distortion  associated  with  the  mapping  from  the  ellip¬ 
soid  to  the  conformal  sphere  over  geographical  areas  that  are 
latltudinally  centered  at  0  and  40  degrees  respectively. 

A  common  feature  of  Figures  6  through  8  Is  that  the  magnifica¬ 
tion  is  greater  than  or  equal  to  unity  at  the  projection  center  and 
increases  as  the  distance  from  the  projection  center  increases. 

What  would  appear  more  desirable,  as  it  would  tend  to  decrease  the 
registration  error  between  observations  of  the  same  aircraft  by 
different  radars,  is  a  unity  magnification  that  would  occur  approxi¬ 
mately  midway  between  the  projection  center  and  the  extremity  of 
radar  coverage.  This  would  effect  magnifications  of  less  than 
unity  at  the  projection  center,  a  contour  of  unity  magnification 
within  the  geographic  area  of  coverage  and  a  magnification  contour 
that  exceeds  unity  at  the  extremity  of  radar  coverage. 

To  accomplish  the  desired  result  equation  (14) ,  which  for  con¬ 
venience  is  reiterated  below,  will  be  considered  when  setting  the 
magnification  k  to  unity  at  a  particular  value  of  (L,X)  and  for  a 
particular  projection  center  (Lo,Xq). 


.  a  2E  Cos  $ 

t  "  N  COS  L  [1  +  Sin  4>  Sin  $  +  Cos  *  Cos  <t>  Cos  (X-X  )] 

o  o  o 
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Solving  equation  (IS)  for  the  resultant  radius  of  the  confornal  sphere 
E,  ve  obtain: 

N  Cos  L  [1  +  Sin  $  Sin  $  +  Cos  $  Cos  $  Cos  (X-X  ) 
rj  „  _ o _ o  o 

2  Cos  $ 


Equation  (16)  determines  the  radius  of  the  conformal  sphere  that  will 
result  in  unity  magnification  at  the  point  (L,X)  for  the  particular 
projection  center  (L  ,X  ).  Table  III  snows  the  radii  of  conformal 
spheres  that  result  from  solving  equation  (16)  for  several  values 
of  (L,X)  when  the  projection  center  (L  ,X  )  is  chosen  as  (40s,  -3s). 
The  initial  value  of  (L,X),  is  of  no  practical  importance  since  the 
projection  center  is  approximately  2400  miles  away  from  the  point 
(0s,  -7s).  Nevertheless,  it  has  been  included  in  Table  III  in  order 
to  better  understand  the  implications  of  equation  (16). 

Table  III 

Radius  of  Conformal  Sphere 

(Ellipsoid  To  Plane) 

L,  X 

Earth  Radius 

i 

o 

3041.62 

42,-7 

3435.45 

38,  0 

3437.21 

40,-3 

3439.29 

Comparing  the  earth  radii  in  Table  III  to  those  in  Table  I,  ve 
note  that  the  earth  radius  corresponding  to  (0,-7)  is  approximately 
400  miles  less  than  any  of  the  other  earth  radii,  tfhen  the  points 
(L,X)  are  proximate  to  the  projection  center,  the  resultant  earth 
radii  are  approximately  equal  to  the  earth  radius  that  was  presented 
in  Table  I  and  which  was  optimized  for  the  mapping  from  the  ellipse 
to  the  sphere.  When  the  point  (L,X)  equals  the  coordinates  of  the 
projection  center,  the  resultant  earth  radius  equals  the  earth  radius 
of  Table  I  which  was  optimized  for  the  mapping  from  the  ellipse  to 
the  sphere , 

♦ 

Figure  9  shows  the  contours  of  equal  magnification  that  result 
from  solving  equation  (15)  while  requiring  that  a  unity  magnification 
be  obtained  at  the  point  0”,  -7s.  The  appropriate  earth  radius  that 
was  used  in  obtaining  th~  contours  is  3041.62  miles  and  the  projection 
center  is  40s,  -3s.  The  geographical  area  over  which  the  contours 
are  presented  is  identical  to  that  of  the  previous  figures. 
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MAGNIFICATION;  L,X»0#,-7* 


Examination  of  Figure  9  shows  that  the  scale  equals  0.99437  at 
the  projection  center  and  increases  with  the  distance  from  the  center. 
If  the  contours  were  plotted  over  an  extended  geographical  area,  the 
magnification  would  eventually  become  unity  at  the  point  0°,  -7°. 

Figure  10  presents  contours  of  equal  magnification  when  L,X 
equals  42°,  -7°  in  equation  (16).  Figure  11  presents  similar  contours 
when  L,X  equals  38°,  0° ,  Finally,  Figure  12  presents  magnification 
contours  for  L,X  equal  to  the  projection  center,  namely  L,X  equals 
40°,  -3°.  Figure  12  is  identical  to  Figure  6. 

Figures  9  through  12  show  that  the  magnification  contours  remain 
symmetric  about  the  vertical  axis  through  the  specified  projection 
center.  Therefore,  to  minimize  the  distortion  associated  with  the 
mapping  from  the  conformal  sphere  to  the  common  coordinate  plane, 
the  projection  center  should  be  located  proximate  to  the  center  of 
the  geographical  area  of  interest. 

Figures  9  through  12  also  show  that  a  unity  magnification  contour 
may  be  obtained  through  a  designated  location  within  the  geographical 
area  of  interest  by  simply  solving  equation  (16)  for  the  appropriate 
radius  of  the  conformal  sphere.  Using  the  resultant  radius  as  the 
value  of  E  in  the  transformation  equations  will  effect  the  desired 
result.  It  is  suggested  that  the  unity  magnification  contour  should 
be  located  approximately  midway  between  the  projection  center  and 
the  maximum  extent  of  radar  coverage  in  order  to  decrease  the  registra¬ 
tion  error  between  observations  of  the  same  aircraft  by  different 
radars . 

EARTH  RADIUS  AND  COORDINATE  CONVERSION 

Having  obtained  the  radius  of  the  conformal  sphere  that  results 
in  a  unity  magnification  along  a  particular  elliptical  contour.  It 
is  necessary  to  determine  whether  this  radius  is  suitable  for  use  in 
the  stereographic  coordinate  conversion  algorithms.  Intuitively  one 
could  argue  that  the  earth's  radius  to  be  used  in  the  coordinate 
conversion  algorithm  should  be  obtained  from  equation  (16)  based  upon 
the  coordinates  of  the  radar  site  location  (L,X)  and  the  origin  of  the 
common  coordinate  plane  (L  ,X  ).  Such  an  approach  would  result  in  a 
unique  earth  radius  being  Sss8ciated  with  each  radar  site. 

To  ascertain  the  effect  of  a  unique  earth  radius  for  each  radar 

site,  the  extent  of  the  variation  in  the  earth's  radius  over  the 

geographical  area  described  by  36"  <  L  £  44*  and  -9°  <_  X  3*  waa 

determined.  Equation  (16)  was  solved  for  varying  (L,X)  within  the 

geographical  area  with  the  origin  of  the  comnon  coordinate  plane 

located  at  L  ■  40"  and  X  -  -3*.  The  resultant  minimum  and  maximum 
o  o 


44n 


values  of  the  earth  radii  were  approximately  3428  and  3440  miles 
respectively.  To  determine  the  impact  of  such  a  variation  in  the 
earth's  radius,  the  stereographic  coordinate  conversion  algorithm 
and  the  impact  of  the  earth's  radius  therein  will  be  considered. 

Stereographic  coordinate  conversion  consists  of  using  slant 
range  and  azimuth  data,  the  height  of  an  aircraft  above  sea  level 
and  the  elevation  of  the  radar  above  sea  level  to  obtain  the  atereo- 
graphic  projection  of  an  aircraft's  location  onto  the  radar  coordinate 
plane.  From  page  31  of  Reference  3,  a  coordinate  conversion  equation 
which  combines  suitable  accuracy  with  minimal  processing  requirements 
is: 


where: 

R  is  the  stereographic  ground  range  in  the  radar  plane. 

S  is  the  measured  slant  range. 

H  is  the  measured  height  of  the  aircraft  above  sea  level. 

h  is  the  elevation  of  the  site  above  sea  level. 

E  Is  the  radius  of  the  earth. 

To  determine  the  variation  In  the  ground  range  R  that  results 
from  using  site  dependent  earth  radii,  equation  (17)  has  been  evaluated 
for  several  values  of  the  numerator  while  E  assumes  the  minimum  and 
maximum  values  for  our  hypothetical  area  of  coverage.  Table  IV  shows 
the  impact  of  using  the  minimum  and  maximum  earth  radii  on  the  com¬ 
putation  of  ground  range. 

Table  IV 


Ground  Range  Versus  Earth  Radius 
R  based  on: 


225  224.8688  224.8693 


100 


99.9417 


99.9419 


Examination  of  Table  IV  show*  that  there  la  no  eeaentlal  differ¬ 
ence  in  the  calculated  ground  range  due  to  the  varying  earth  radii. 
Therefore,  the  earth  radius  as  derived  froa  equation  (16)  and  which 
yields  a  unity  magnification  contour  at  a  location  approximately 
midway  between  the  projection  center  and  the  maximum  extent  of  radar 
coverage  should  also  be  used  in  the  coordinate  conversion  algorithm. 
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SECTION  IV 


SUMMARY  OP  STEREOGRAPHIC  PROJECTION  ALGORITHM 


INTRODUCTION 

This  section  of  the  report  reiterates  the  more  important  results 
that  have  been  indicated  herein.  Additionally,  a  summary  of  the 
stereographic  conversion  and  transformation  equations  is  presented. 
The  conversion  algorithms  have  been  developed  in  Reference  3  and  the 
transformation  equations  have  been  derived  herein. 


ORIGIN  OF  COMMON  COORDINATE  PLANE 

The  latitude  and  longitude  of  the  origin  of  the  common  coordi¬ 
nate  plane  are  stored  parameters  that  are  required  for  stereographic 
projection  of  radar  data.  This  origin  should  be  chosen  proximate  to 
the  center  of  the  geographical  area  of  Interest  to  equalize  the  map¬ 
ping  distortion  about  this  location  and  to  minimize  the  overall 
distortion. 


CHOICE  OP  EARTH'S  RADIUS 

Stereographic  projection  requires  the  mapping  of  points  on  or 
above  a  designated  location  on  the  earth  ellipsoid  onto  a  sphere 
which  is  called  the  conformal  sphere.  The  radius  of  this  sphere  is 
a  stored  parameter  that  is  required  by  joth  the  stereographic  trans¬ 
formation  and  conversion  algorithms. 

Having  selected  the  origin  of  the  common  coordinate  plane  as 
described  above  and  having  selected  a  particular  latitude  and  longi¬ 
tude  through  which  a  unity  magnification  contour  is  desired,  the 
radius  of  the  conformal  sphere  should  be  calculated  from  equation  (18) . 
It  is  recommended  that  the  unity  magnification  contour  should  be 
located  approximately  midway  between  the  origin  of  the  common  coordi¬ 
nate  plane  and  the  maximal  extent  of  the  radar  coverage  in  order  to 
decrease  the  registration  error  between  observations  of  the  same  air¬ 
craft  by  different  radars.  The  radius  of  the  conformal  sphere  that 
is  obtained  from  equation  (18)  should  be  used  in  both  the  stereo¬ 
graphic  conversion  and  transformation  algorithms. 
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(18) 


N  Cos  L  [1  +  Sin  $  Sin  $  +  Cos  $  Cos  <p  Cos  (X-X  )] 

o  o  o 

2  Cos  $ 


where: 

E  is  the  radius  of  the  conformal  sphere. 

N  is  Eq(l-e2Sin2L)’Js  where  Eq  -  3444.054  miles  Is 

the  earth' 8  equatorial  radius  and  where  e2  ■  .00672267 
is  the  earth's  eccentricity  squared. 

L  is  the  geographic  latitude  of  a  designated  loca¬ 

tion.  This  is  the  angle  which  a  normal  to  the 
surface  of  an  ellipsoid  makes  with  the  equatorial 
plane  as  shown  in  Figure  4. 

4>,X  is  the  conformal  latitude  and  longitude  of  a 

point  through  which  the  unity  magnification  con¬ 
tour  is  desired.  $  is  obtained  via  equation  (19). 

<j>o,XQ  is  the  conformal  latitude  and  longitude  of  the 
origin  of  the  common  coordinate  plane.  4>0  is 
obtained  via  equation  (19) . 

CONFORMAL  AND  GEOGRAPHIC  LATITUDE 

The  conformal  latitude  may  be  obtained  from  the  geographic 
latitude  according  to  the  relationship: 


Tan 


Tan 


[(*)  ♦  (*)]  ’  [ 


1  -  e  Sin  L 
1  +  e  Sin  L 


(19) 


where: 

<fr  is  the  conformal  latitude  of  a  point,  the  geographic 
latitude  of  which  is  L. 

e  •  .08199189  is  the  earth's  eccentricity. 

For  computational  convenience,  the  conformal  latitude  $  may  be  obtained 
from  the  geographic  latitude  L  by  a  series  approximation  to  equation  (19) . 


STEREOGRAPHIC  COORDINATE  CONVERSION 


Slant  range,  azimuth  and  height  data  should  be  projected  Into 
the  radar  tangent  plane  as  described  in  Reference  3.  When  height 
data  is  available  the  stereographic  ground  range  R,  the  distance 
from  the  radar  to  the  aircraft  in  the  radar  plane,  should  be  cal¬ 
culated  as  follows: 


.  (S2  -  tt-w2]1* 

1  +  (2/E) 

When  height  data  is  unavailable,  the  stereographic  ground  range  R 
should  be  calculated  according  to: 


R  -  1.0025  S  -  0.65 

where: 

R  is  the  stereographic  ground  range. 

S  is  the  measured  slant  range. 

H  is  the  measured  aircraft  altitude  above  sea  level, 
h  is  the  radar  site  elevation  above  sea  level. 

E  is  the  radius  of  the  conformal  sphere  from  equation  (18) . 


(21) 


RADAR  SITE  COORDINATES 

The  rectangular  coordinates  of  a  radar  site's  location  or  any 
other  point  on  the  earth's  surface  relative  to  the  origin  of  the  com¬ 
mon  coordinate  plane  should  be  calculated  from  equations  (22) .  Posi¬ 
tive  values  of  Ur  and  Vr  correspond  to  locations  that  are  east  and 
north  of  the  origin  of  the  connon  coordinate  plane. 
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(22) 


Sin  (X-A  )  Cos  $ 

2E  ° 

1  +  Sin  <j)  Sin  4  +  Cos  $  Cos  $  Cos  (X-X  ) 

°  o  o 

Sin  <j>  Cos  <i  -  Cos  <J>  Sin  <j»  Cos  (X-X  ) 

2E  ^ - - - 2 _ o 

1  +  Sin  $  Sin  4>q  +  Cos  $  Cos  <j>o  Cos  (A-Xq) 

ur»  vr  are  the  rectangular  coordinates  of  the  point  designated 
by  (L,X)  in  the  common  coordinate  plane. 

‘J1 » ^  are  the  conformal  latitude  and  longitude  of  the  point 

to  be  projected. 

are  the  conformal  latitude  and  longitude  of  the  origin 
of  the  common  coordinate  plane. 

E  ia  the  radius  of  the  conformal  sphere  fro*  equation  (18) . 

Latitudes  north  of  the  equator  are  defined  as  positive  whereas 
latitudes  south  of  the  equator  are  negative.  Longitudes  eacs.  oi  the 
prime  meridian  are  defined  as  positive  and  those  west  of  the  prime 
meridian  are  negative. 


U 


where : 


ANGULAR  ADJUSTMENT  OF  RADAR  DATA 

Radar  data,  which  are  measured  relative  to  true  north  at  the 
radar  site  location,  should  be  adjusted  so  that  the  azimuth  is  rela¬ 
tive  to  north  at  the  origin  of  the  common  coordinate  system.  The 
angular  adjustment  8  is  defined  as  follows: 

(Sin  0  +  Sin  $o)  Sin  (X-X) 

Cos  *  Cos  <|>o  +  (1  +  Sin  *  Sin  4>q)  Cos  (X-Xq) 


(23) 
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where; 


4> * X  are  the  conformal  latitude  and  longitude  of  the 

location  of  the  radar  site, 

$  ,A  are  the  conformal  latitude  and  longitude  of  the 
origin  of  the  common  coordinate  plane. 

An  azimuth  datum  6,  which  is  measured  relative  to  true  north  at 
the  radar  site,  should  be  adjusted  to  true  north  at  the  origin  of  the 
common  coordinate  system  as  follows: 


9  -6  +  9 


TRANSFORMATION  FROM  RADAR  TO  COMMON  COORDINATES 

The  transformation  from  the  radar  coordinate  plane  to  the  common 
coordinate  plane  should  be  accomplished  by  an  approximation  to  the 
stereographic  transformation  equation. 

When  aircraft  height  information  is  available,  the  rectangular 
coordinates  of  an  aircraft  with  respect  to  the  common  coordinate 
origin  should  be  calculated  as  follows: 


U  -  U  +  K 
r 


R  Sin  (0+0)  + 


RZ|Wr|Sin  [2(0+6)  -  y] 


4E 


V  +  K 
r 


R  Cos  (8+6)  + 


R  W  | Cos  [2(0+6)  -  y) 


4E 


where: 

Ur,Vr  are  the  rectangular  coordinates  of  the  radar 
site  in  the  cosnon  coordinate  plane. 

R  is  the  stereographic  ground  range  in  the  radar 

plane  from  equation  (20) , 


K 


-  1  +  |Wr|2/4E2 

l«rl  -  iur2  +  v.2]*1 

Y  “  Tan"1  (TTr/Vr) 

Equations  (24)  nay  be  calculated  by  the  following  expressions  which 
preserve  cooputational  accuracy  while  ninlntsing  the  processing 
requirements . 


A2|w  I  Sin  [2(0+6)  -  y! 

U  *  U  +  T 

A  Sin  (0+6)  +  — ^ - n - 

r 

4E2 

A2|W  1  Cos  [2(9+6)  -  y! 

V  »  V  +  T 

A  Cos  (6+6)  +  — - 5 - 

r 

4EZ 

where : 


T 


2  2  h 

A  -  ISZ  -  (H-h;  j 

When  aircraft  height  infomation  is  unavailable,  the  rectangular 
coordinates  of  an  aircraft  with  respect  to  the  common  coordinate 
origin  should  be  calculated  from; 


U  -  Ur  +  K  R  Sin  (9+6) 

V  -  V  +  K  R  Cos  (0+6) 
r 


(26) 


where: 

R  is  the  stereographic  ground  range  In  the  radar  plane  fron 
equation  (21). 
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APPENDIX  A 


THE  STEREOGRAPRIC  TRANSFORMATION  EQUATION 


Fron  pages  86  and  133  of  Reference  5,  a  point  on  the  surface  of 
a  spherical  earth,  the  latitude  of  which  is  L  and  the  longitude  of 
which  is  X,  is  stereo graphically  projected  onto  a  plane  with  an  origin 
at  a  latitude  of  LQ  and  a  longitude  of  XQ  by  equation  (A-l) .  The 
resultant  point  in  the  plane  of  projection  la  designated  by  Z  and  its 
rectangular  components  are  X  and  Y. 


Z  -  X  +  iY 


exp  [*s(t  -iX-t o+1Xq)  1  -  exp  [-^(T-iX-Tp+iX^J 

exp  [*i  (t-IX+t  +i  X  )  ]  +  exp  [-^(t-IXHK  +±X  )] 
r  0  0  o  o 


where: 


eT  -  Tan  (ir/4  +  L/2) 

T 

e  °  »  Tan  (ir/4  +  Lq/2) 

L  “  latitude  of  the  point  to  be  projected 

Lq  •  latitude  of  the  coordinate  origin 

X  *  longitude  of  the  point  to  be  projected 

Xq  «  longitude  of  the  coordinate  origin 

E  -  radius  of  a  spherical  earth 


Substituting  A  -  ^(t-IX),  B  -  *s(-t  +iX  )  and  C  -  %(x  4-iX  ). 

O  o  o  o 

D  m  ^(-T^+iX^)  and  F  -  *s(t^-KX^),  Z  nay  be  expressed  as: 


2E1 


A  B  -A  -B 
e  e  -  e  e 

A  C  ,  -A  -C 
e  e  +  e  e 


-  2Ei 


eB  -  e~2A  e-B 

C  ,  -2A  -C 

e  +  e  e 
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(A-l) 


(A-2) 


Preceding  pigs  blink 


Solving  equation  (A-2)  for  exp  (t-1X) ,  we  obtain 


Z  +  2Ei  exp  (t  ) 

exp  (2A)  *  exp  (r-iX)  -  2Et  exp  '(1X~)'~-  Z  exp  (t  Tl T) 

o  o  o 

Let  us  now  stereographically  project  the  point  (L,X)  onto  a 
plane  with,  an  origin  of  (L^,X^>.  The  resultant  point  in  the  plane 
of  projection  Is  designated  by  V  and  Its  rectangular  coordinates 
are  U  and  V.  The  form  of  W  will  be  identical  to  equation  (A— 1) 
except  that  tq  will  be  replaced  by  and  X0  will  be  replaced  by  X^. 
The  point  (L^,X^)  may  be  considered  as  the  origin  of  the  common 
coordinate  plane  and  the  point  (1^,  ,XQ)  may  be  regarded  as  the  origin 
of  the  radar  plane.  Rewriting  equation  (A-l)  as  suggested,  W  is: 


(A-4) 


Substituting  from  (A-3)  into  (A-4)  and  performing  the  requisite 
manipulations,  V  may  be  expressed  as: 


Z  exp(-x1 )+2Ei  exp (t  -t. )-2Ei  exp[i(X  -X-)]+Z  exp[x  +i(X  -X, ) ] 

2V±  _ i _ 0  -*• _  0  1.  0  O  1 

Z+2E1  exp(xo)+2Ei  exp[-f1+i(Xo-X1)]-Z  exp[  C^-T^Hi^-X^] 

Z{ exp [ i(XQ-X  ] +exp [ - (tq+t  ] }  +2Ei{ exp (-t^) -exp [ -xo+i(Xo-X^) J } 

1  +  «q>[“(To-K1)+l(lo-X1)3  -  Z  2El{exp(-To)-exp[-T1-H(X^-X1)]} 

4E2 

Therefore,  H  may  be  written  as: 


ft , 

exp  [i(Xo-X1>]  +  exp  [-(to+t1)J 

lexpC-,!)  - 

L 

1+exp  [-(x^^  + 

j  1+exp  [-(^p+^j^)  | 

1  ^  “jv* 

exp  ("T0>  ”  expf-Tj+KX^-?^)] 

4E 

1  +  exp  [-(t^+Tj)  +i(XQ-Xl)J 
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Equation  (A-5)  is  the  stereographic  transformation  equation. 

We  will  now  examine  this  equation  term  by  term  in  order  to  express 
it  in  a  simpler  form. 

The  fraction  by  which  Z  is  multiplied  in  the  numerator  of  equa¬ 
tion  (A-5)  may  be  written  as: 


exp  [i(AQ-X1)]  +  exp 
1  +  exp  C-(to-Ht1)  +  i(XQ-X1)] 


(eT  +  Cos  A)  +  i  Sin  A 
(1  +  eT  Cos  A)  +  ieT  Sin  A 


where :  t 

A 


j 


-  (Xo  -  Xl> 

■  Tan  [Js(x/2  +  L^] 

-  Tan  [i*(x/2  -  Lj] 


fl  +  Sin  L,"1 

f  Cos  L.  1 

i 

1 

L  Cos  Lj 

.1  -  Sin  Lj  j 

[l  "  Sln  LJ 

c°’Li  1 

L  Cos 

1  +  Sin  Lj  J 

*  0  or  1 


Multiplying  the  numerator  and  denominator  of  (A-6)  by  It a  com¬ 
plex  conjugate,  and  performing  the  necessary  algebra,  (A-6)  may  be 
expressed  as: 


2eT  +  e2x  Cos  A  +  Cos  A  -f-  i[Sin  A  (l-e2*)] 
1  +  2eT  Cos  A  +  e^T 


(A-6) 


(A-7) 


Calculating  the  real  (R)  and  imaginary  (I)  portions  of  expression 
(A-7)  via  the  above  definitions  of  exp  (t.)  and  exp  (-t.),  R  and  I  may 
be  written  as:  J  3 


Cos  Lq  Cos  Lx  +  (1  +  Sin  Lq  Sin  1^)  Cos  (A^-Aj) 
•  1  +  Sin  Lq  Sin  +  Cos  Lq  Cos  Lj  Cos  (A  -Xj} 

(Sin  Lq  +  3 in  1^)  Sin  (X  -A  ) 

1  +  Sin  Lq  Sin  Lx  +  Cos  Cos  1^  Cos  (X^A^ 


(A-8) 
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2  2 

Although  tedious  to  verify,  it  may  be  shown  that  R  +1  equals 
unity.  Consequently,  (A-6)  may  be  expressed  In  terms  of  its 
polar  coordinates  as  follows: 


exp  [-<tq  +  tl)]  +  exp  [iC^-Xj^)]  =  [r2*!2]5*  g-iB  _  e'iP 

1  +  exp  [-(tq  +  t1)  +  i(Xo-X1)] 

where 

(Sin  Lq  +  Sin  Lj)  Sin  (X-j-A^ 

Cos  L  Cos  L.  +  (1  +  Sin  L  Sin  L1 )  Cos  (X  -A) 

O  1  OX  O  X 

and 

L  ,X  are  the  latitude  and  longitude  of  the  origin  of 
the  radar  coordinate  plane. 

1*1 ’^1  are  the  latitude  and  longitude  of  the  origin  of 
the  common  coordinate  plane. 

From  consideration  of  the  definition  of  B,  it  will  be  noted  that 
a  positive  3  angle  is  measured  in  a  counterclockwise  direction  from 
the  R  axis.  Having  reduced  (A-6)  to  the  leftmost  term 

in  the  numerator  of  equation  (A-5)  may  be  expressed  by  2e-^ . 

The  right  side  of  the  numerator  of  equation  (A-5)  is: 

exp  (— t.)  -  exp  [— t  +  i(X  -X.)J 

JTfJ  - - 1 - - U  - 

1  +  exp  [-frQ  +  t^)  +  i(XQ-X1)] 


exp(To-iXo)  -  exp(T^-iXj.)  .+  exp(-T1~iX1)  -  erp(-2iXj-To+iAo) 
[expC-iXj^)  +  exp(TQ+T1-iXo)]  {l  +  exp  [-(i^-h^)  +  i(XQ-X1)]) 


2E1 


exp  (Tq-IAq)  -  exp  (r^-iX^) 
exp  (-iA^  +  exp  (T0+T1-1X0) 


exp  "  “P  [  Js(T0~iX0“Tl+i*l^ 

exp  [35(To-iXo+r1+iX1)]  -t-  exp  [-^(^-iA^j+iXp] 


W 


r 


(A-9) 


(A-10) 


Therefore,  the  rlgnt  side  of  the  numerator  of  equation  (A-5)  is  W_ , 
which  represents  the  origin  of  the  radar  plane  projected  onto  the 
common  coordinate  plane.  The  entire  numerator  of  equation  (A-5)  is 
(Ze-i8+wr). 

The  fraction  by  which  (-Z/4E2)  is  multiplied  in  the  denominator 
of  equation  (A-5)  is; 


exp  (— T  )  -  exp  [-t  ,+i(X  -A  )]  exp  (-T  )  -  exp  [-T.+i(X  -X.)] 

2Ei  - - - - - 2_i —  -  2Ei  - - - ± _ ± 

1  +  exp  l-(VTl>  +  i<VXl)1  £Xp  UOq-V1  +  “*»  I-(To+Tl)1 

exp  [i(XQ-X1)]  +  exp  [-(t^-Kj)] 

1  +  exp  f-<To+r1)  +  i(XQ-X1)] 


which  from  expression  (A-9)  is: 


2Ei 


exp  (-to)  -  exp  [-Tj  +  i(XQ-X1)] 


-10 


exp  (i(Xo-X1)]  H  exp  [-(t^+Tj)] 

The  terra  within  the  parenthesis  may  be  expressed  as: 

exp  (xo+i(2Xo-X1)]  +  exp  (-Tj+iX^  -  exp  (tj+iAo)  -  exp  (-t^iXj) 
— 2Ei  " 

[exp  (to+t1+1Xo)  +  exp  (iX^]  {  exp  [1(Xq-X1)]  +  exp  [-(t^+Tj)]} 


(A-ll) 


-2Ei 


exp  (tq  +  iXQ)  -  exp  (tl  +  iAj) 
exp  (t  +  t.  +  IX  )  +  exp  (IX.) 


-2Ei 


exp  {«j(To+i>o-T1-iX1)3  -  exp  [-M^+iX^^-iX^] 
exp  tJS^t0+lx0+^1”iix>3  ♦  exP  t Js^T0+iXo'Kl_iXl^ 


-  W 


(AnL2) 
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where  W  -  Ur  -  iVr  Is  the  complex  conjugate  of  Wr.  Therefore,  ’.he 
fractioi  by  which  (-Z/4E*)  is  multiplied  in  the  denominator  of  equa¬ 
tion  (A-5)  is  W  e”*®. 

r 

Substituting  the  results  of  expressions  (A-9)  through  (A-12)  into 
equation  (A-5)  yields! 


W 


Ze"1S  +  W 

-TTT-fe 

r 


4E 


(A-13) 


APPENDIX  B 


PROJECTING  THE  EARTH'S  SURFACE  ONTO  A  PLANE 


This  Appendix  derives  for  a  spherical  earth  the  equations  which 
stereographically  project  points  on  the  surface  of  the  earth  onto  a 
plane  of  rectangular  coordinates.  Figure  B-l  depicts  the  geometry 
for  a  spherical  triangle  with  the  origin  of  the  coordinates  at  lati¬ 
tude  L  and  longitude  XQ.  The  ooint  to  be  projected  has  coordinates 
L.X . 


In  Figure  B-l,  N  and  E  are  the  directions  of  north  and  east 
respectively  in  the  plane  of  projection.  The  line  OP'  is  the  arc 
OF  projected  into  the  plane.  The  angle  y  in  the  plane  corresponds 
to  the  spherical  angle  ZOP .  0  is  the  angle  between  the  radii  drawn 

from  the  earth's  center  to  the  origin  of  the  coordinate  system  and 
to  the  point  to  be  projected. 

From  the  law  of  cosines  for  the  spherical  triangle  OPZ,  we  obtain: 

Cos  PO  -  Cos  PZ  Cos  OZ  +  Sin  PZ  Sin  OZ  Cos  PZO 

(B-l) 

Cos  0  -  Sin  L  Sin  L  +  Cos  L  Cos  L  Cos  (X-X  ) 

o  o  o 

From  spherical  triginometry,  the  relationship  between  two  angles 
and  three  sides  is: 

Sin  PO  Cos  ZOP  ■  Cos  PZ  Sin  ZO  -  Cos  ZO  Sin  PZ  Cos  PZO 

(B— 2  ) 

Sin  0  Cos  y  ■  Sin  L  Cos  Lq  -  Sin  Lq  Cos  L  Cos  (X-Xq) 

From  the  law  of  sines : 

Sin  PO  Sin  PZ  Sin  0  Cos  L 

Sin  PZO  ”  Sin  ZOP  “  Sin  (X-X  )  "  Sin  y 

o 

Therefore,  Sin  0  Sin  y  ■  Sin  (X-X  )  Cos  L  (B-3) 

From  equation  (2)  of  Reference  3,  the  stereographic  ground  range 
designated  |Wr|  herein,  may  be  obtained  from  the  relationship: 

IwJ  -  2E  Tan  (0/2)  -  2E  i  ■, 
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EARTH 

CENTER 


The  ur»Vr  components  of  the  projected  point  Wr  are  obtained  from: 


U 


r 


Wr|  Sin  y 


2E  Sin  0  Sin  y 
1  +  Cos  0 


V 


r 


Cos  y 


2E  Sin  0  Cos  y 
1  +  Cos  0 


(B-4) 


Substituting  from  equations  (B-l) ,  (B-2) ,  and  (B-3)  into  (B-4) ,  we 
obtain: 

2E  Sin  (X-X  )  Cos  L 

n  .  .  . ° . - . 

r  1  +  Sin  h  Sin  L  +  Cos  L  Cos  L  Cos  (X-X  ) 

o  op 


2E[Sin  L  Cos  L  -  Sin  L  Cos  L  Cos  (X-XJ] 
o  o  o 

Vr  =  1  +  Sin  L  Sin  L  +  Cos  L  Cos  L  Cos  (X-X) 

o  o  o 

where : 

L , X  are  the  latitude  and  longitude  of  the  point  to 

be  projected  onto  the  coordinate  plane. 

L  , X  are  the  latitude  and  longitude  of  the  origin 

°  of  the  plane  of  rectangular  coordinates. 

U  , V  are  the  rectangular  coordinates  of  the  point 

r  L,X  in  the  plane  of  projection. 


(B-5) 
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